Application of factor analytic and spatial mixed models for the analysis of multi-environment trials in common bean (Phaseolus vulgaris L.) in Ethiopia.

Common bean (Phaseolus vulgaris L.) is one of the most important grain legumes consumed globally, especially in Ethiopia, for its edible seeds, cash crops, and supply of protein for farmers. Efficient statistical methods must be employed for the evaluation of common bean varieties to accurately select superior varieties that contribute to agricultural productivity. The objective of this study was to identify promising large mottled bean varieties through analysis of multi-environment trials (MET) data using multiplicative spatial mixed models. In this study, 16-18 large mottled common bean varieties, including one check, were sown across nine growing environments in Ethiopia using lattice and alpha lattice designs, with three replications laid out in a square or rectangular (row by column) array of plots, respectively during the main cropping season from 2015 to 2018. We present a linear mixed model analysis that integrates spatial and factor analytic (FA) models, and the heritability measure was used to evaluate the efficiency of these models with the conventional analysis. The analysis of the spatial model, and more significantly, the spatial+FA model, revealed a notable enhancement in heritability. With the exception of a trial conducted at Kobo, a genotype DAP 292, found to be good performing for days to flowering and maturity, but for yield only across four clusters of trials, C2, C3, C5 and C7, formed with trials of relatively high genetic variance. Across these clusters, the yield advantage of this variety over the check ranged from 10-32%. This genotype also has a yield that is somewhat comparable to the check across the remaining clusters. Overall, both the spatial and factor analytic models proved to be effective approaches for analyzing the data in this study. The analysis of multi-environment trial data through the use of more efficient statistical models can provide a more robust platform for evaluating common bean varieties with greater confidence in selecting superior varieties across a range of environments. Hence, scaling up the use of this efficient analysis method is indispensable for enhancing the selection of superior varieties.

Phenotypic and yield responses of common bean (Phaseolus vulgaris l.) varieties to different soil moisture levels.

BACKGROUND: Morphological plasticity is one of the capacities of plants to modify their morphological appearance in response to external stimuli. A plant's morphology and physiology are constantly tuned to its variable surroundings by complex interactions between environmental stimuli and internal signals. In most of plant species,, such phenotypic and physiological expression varies among different varieties based on their levels of particular environmental stress conditions. However, the morphological and yield responses of common bean varieties to different environmental conditions are not well known. The purpose of the study was to evaluate morphological and yield response of common bean to soil moisture stress and to investigate the morphological mechanism by which common bean varieties tolerate fluctuations in moisture stress. METHODS: A pot experiment was carried out to investigate the effects of different moisture levels on the phenotypic and yield responses of common bean varieties. A factorial combination of five common bean varieties (Hirna, kufanzik, Awash-1, Ado, and Chercher) and three moisture levels (control, waterlogging stress, and moisture deficit stress) was used in three replications. Moisture stress treatments were started 20 days after planting, at the trifoliate growth stage. To evaluate the response of each variety, morphological and yield data were collected at week intervals. MAIN RESULTS: The results indicated that moisture levels and varieties had a significant influence on all growth parameters. Crop phenology was significantly influenced by the interaction effect of moisture level and variety. Exposing Hirna variety to moisture stress led to extended flowering and pod setting by 23 and 24 days, respectively, compared to the other treatments. The results showed that the phenotypic responses to moisture deficit and waterlogging stress varied between varieties. Waterlogging stress had a stronger reduction effect on the fresh weight, dry weight and leaf area of common bean varieties than moisture deficit and the control. Pods per plant, seeds per plant, grain yield per plant, and harvest index were significantly influenced by the varieties, moisture stress levels and their interaction. Except for Chercher and Hirna. However, varieties Ado, kufanzik and Awasha-1 did not show significant differences on the time of flower initiation due to moisture level. Biomass and growth in leaf fresh weight, leaf dry weight, leaf area, leaf number and plant height were significantly influenced by moisture level. When moisture deficit and waterlogging stress occurred, Ado and Awash-1 were more responsive to moisture stress than Hirna, Chercher, and Kufanzik. CONCLUSION: Hence, Hirna and Kufanzik varieties were found to be tolerant because they produced higher yields than the Chercher, Awash-1, and Ado varieties.

A Past Genetic Bottleneck from Argentine Beans and a Selective Sweep Led to the Race Chile of the Common Bean (Phaseolus vulgaris L.).

The domestication process of the common bean gave rise to six different races which come from the two ancestral genetic pools, the Mesoamerican (Durango, Jalisco, and Mesoamerica races) and the Andean (New Granada, Peru, and Chile races). In this study, a collection of 281 common bean landraces from Chile was analyzed using a 12K-SNP microarray. Additionally, 401 accessions representing the rest of the five common bean races were analyzed. A total of 2543 SNPs allowed us to differentiate a genetic group of 165 accessions that corresponds to the race Chile, 90 of which were classified as pure accessions, such as the bean types 'Tórtola', 'Sapito', 'Coscorrón', and 'Frutilla'. Our genetic analysis indicates that the race Chile has a close relationship with accessions from Argentina, suggesting that nomadic ancestral peoples introduced the bean seed to Chile. Previous archaeological and genetic studies support this hypothesis. Additionally, the low genetic diversity (π = 0.053; uHe = 0.53) and the negative value of Tajima' D (D = -1.371) indicate that the race Chile suffered a bottleneck and a selective sweep after its introduction, supporting the hypothesis that a small group of Argentine bean genotypes led to the race Chile. A total of 235 genes were identified within haplotype blocks detected exclusively in the race Chile, most of them involved in signal transduction, supporting the hypothesis that intracellular signaling pathways play a fundamental role in the adaptation of organisms to changes in the environment. To date, our findings are the most complete investigation associated with the origin of the race Chile of common bean.

GATA transcription factor in common bean: A comprehensive genome-wide functional characterization, identification, and abiotic stress response evaluation.

The GATA transcription factors (TFs) have been extensively studied for its regulatory role in various biological processes in many plant species. The functional and molecular mechanism of GATA TFs in regulating tolerance to abiotic stress has not yet been studied in the common bean. This study analyzed the functional identity of the GATA gene family in the P. vulgaris genome under different abiotic and phytohormonal stress. The GATA gene family was systematically investigated in the P. vulgaris genome, and 31 PvGATA TFs were identified. The study found that 18 out of 31 PvGATA genes had undergone duplication events, emphasizing the role of gene duplication in GATA gene expansion. All the PvGATA genes were classified into four significant subfamilies, with 8, 3, 6, and 13 members in each subfamily (subfamilies I, II, III, and IV), respectively. All PvGATA protein sequences contained a single GATA domain, but subfamily II members had additional domains such as CCT and tify. A total of 799 promoter cis-regulatory elements (CREs) were predicted in the PvGATAs. Additionally, we used qRT-PCR to investigate the expression profiles of five PvGATA genes in the common bean roots under abiotic conditions. The results suggest that PvGATA01/10/25/28 may play crucial roles in regulating plant resistance against salt and drought stress and may be involved in phytohormone-mediated stress signaling pathways. PvGATA28 was selected for overexpression and cloned into N. benthamiana using Agrobacterium-mediated transformation. Transgenic lines were subjected to abiotic stress, and results showed a significant tolerance of transgenic lines to stress conditions compared to wild-type counterparts. The seed germination assay suggested an extended dormancy of transgenic lines compared to wild-type lines. This study provides a comprehensive analysis of the PvGATA gene family, which can serve as a foundation for future research on the function of GATA TFs in abiotic stress tolerance in common bean plants.

Investigating the genetic basis of salt-tolerance in common bean: a genome-wide association study at the early vegetative stage.

Salinity poses a significant challenge to global crop productivity, affecting approximately 20% of cultivated and 33% of irrigated farmland, and this issue is on the rise. Negative impact of salinity on plant development and metabolism leads to physiological and morphological alterations mainly due to high ion concentration in tissues and the reduced water and nutrients uptake. Common bean (Phaseolus vulgaris L.), a staple food crop accounting for a substantial portion of consumed grain legumes worldwide, is highly susceptible to salt stress resulting in noticeable reduction in dry matter gain in roots and shoots even at low salt concentrations. In this study we screened a common bean panel of diversity encompassing 192 homozygous genotypes for salt tolerance at seedling stage. Phenotypic data were leveraged to identify genomic regions involved in salt stress tolerance in the species through GWAS. We detected seven significant associations between shoot dry weight and SNP markers. The candidate genes, in linkage with the regions associated to salt tolerance or harbouring the detected SNP, showed strong homology with genes known to be involved in salt tolerance in Arabidopsis. Our findings provide valuable insights onto the genetic control of salt tolerance in common bean and represent a first contribution to address the challenge of salinity-induced yield losses in this species and poses the ground to eventually breed salt tolerant common bean varieties.

Genome-wide identification and comparative analysis of the Amino Acid Transporter (AAT) gene family and their roles during Phaseolus vulgaris symbioses.

Amino acid transporters (AATs) are essential integral membrane proteins that serve multiple roles, such as facilitating the transport of amino acids across cell membranes. They play a crucial role in the growth and development of plants. Phaseolus vulgaris, a significant legume crop, serves as a valuable model for studying root symbiosis. In this study, we have conducted an exploration of the AAT gene family in P. vulgaris. In this research, we identified 84 AAT genes within the P. vulgaris genome sequence and categorized them into 12 subfamilies based on their similarity and phylogenetic relationships with AATs found in Arabidopsis and rice. Interestingly, these AAT genes were not evenly distributed across the chromosomes of P. vulgaris . Instead, there was an unusual concentration of these genes located toward the outer edges of chromosomal arms. Upon conducting motif analysis and gene structural analysis, we observed a consistent presence of similar motifs and an intron-exon distribution pattern among the subfamilies. When we analyzed the expression profiles of PvAAT genes, we noted tissue-specific expression patterns. Furthermore, our investigation into AAT gene expression under rhizobial and mycorrhizal symbiotic conditions revealed that certain genes exhibited high levels of expression. Specifically, ATLa5 and LHT2 was notably upregulated under both symbiotic conditions. These findings point towards a potential role of AATs in the context of rhizobial and mycorrhizal symbiosis in P. vulgaris, in addition to their well-established regulatory functions.

Simulations of multiple breeding strategy scenarios in common bean for assessing genomic selection accuracy and model updating.

The aim of this study was to evaluate the accuracy of the ridge regression best linear unbiased prediction model across different traits, parent population sizes, and breeding strategies when estimating breeding values in common bean (Phaseolus vulgaris). Genomic selection was implemented to make selections within a breeding cycle and compared across five different breeding strategies (single seed descent, mass selection, pedigree method, modified pedigree method, and bulk breeding) following 10 breeding cycles. The model was trained on a simulated population of recombinant inbreds genotyped for 1010 single nucleotide polymorphism markers including 38 known quantitative trait loci identified in the literature. These QTL included 11 for seed yield, eight for white mold disease incidence, and 19 for days to flowering. Simulation results revealed that realized accuracies fluctuate depending on the factors investigated: trait genetic architecture, breeding strategy, and the number of initial parents used to begin the first breeding cycle. Trait architecture and breeding strategy appeared to have a larger impact on accuracy than the initial number of parents. Generally, maximum accuracies (in terms of the correlation between true and estimated breeding value) were consistently achieved under a mass selection strategy, pedigree method, and single seed descent method depending on the simulation parameters being tested. This study also investigated model updating, which involves retraining the prediction model with a new set of genotypes and phenotypes that have a closer relation to the population being tested. While it has been repeatedly shown that model updating generally improves prediction accuracy, it benefited some breeding strategies more than others. For low heritability traits (e.g., yield), conventional phenotype-based selection methods showed consistent rates of genetic gain, but genetic gain under genomic selection reached a plateau after fewer cycles. This plateauing is likely a cause of faster fixation of alleles and a diminishing of genetic variance when selections are made based on estimated breeding value as opposed to phenotype.

A comprehensive, improved protocol for generating common bean (Phaseolus vulgaris L.) transgenic hairy roots and their use in reverse-genetics studies.

Generating transgenic hairy roots has been the preferred strategy for molecular studies in common bean (Phaseolus vulgaris L.), since generating stable knockout lines in this species is challenging. However, the number of plants producing hairy roots following the original protocol published in 2007 is usually low, which has impeded progress. Since its initial publication, the original protocol has been extensively modified, but these modifications have not been adequately or systematically reported, making it difficult to assess the reproducibility of the method. The protocol presented here is an update and expansion of the original method. Importantly, it includes new, critical steps for generating transgenic hairy roots and using them in molecular analyses based on reverse-genetics approaches. Using this protocol, the expression of two different genes, used as an example, was significantly increased or decreased in approximately 30% of the transformed plants. In addition, the promoter activity of a given gene was observed, and the infection process of rhizobia in transgenic hairy roots was monitored successfully. Thus, this improved protocol can be used to upregulate, downregulate, and perform promoter activity analysis of various genes in common bean transgenic hairy roots as well as to track rhizobia infection.

Mapping of adult plant recessive resistance to anthracnose in Indian common bean landrace Baspa/KRC 8.

BACKGROUND: The common bean (Phaseolus vulgaris) has become the food of choice owing to its wealthy nutritional profile, leading to a considerable increase in its cultivation worldwide. However, anthracnose has been a major impediment to production and productivity, as elite bean cultivars are vulnerable to this disease. To overcome barriers in crop production, scientists worldwide are working towards enhancing the genetic diversity of crops. One way to achieve this is by introducing novel genes from related crops, including landraces like KRC 8. This particular landrace, found in the North Western Himalayan region, has shown adult plant resistance against anthracnose and also possesses a recessive resistance gene. METHODS AND RESULTS: In this study, a population of 179 F2:9 RIL individuals (Jawala × KRC 8) was evaluated at both phenotypic and genotypic levels using over 830 diverse molecular markers to map the resistance gene present in KRC 8. We have successfully mapped a resistance gene to chromosome Pv01 using four SSR markers, namely IAC 238, IAC 235, IAC 259, and BM 146. The marker IAC 238 is closely linked to the gene with a distance of 0.29 cM, while the other markers flank the recessive resistance gene at 10.87 cM (IAC 259), 17.80 cM (BM 146), and 25.22 cM (IAC 235). Previously, a single recessive anthracnose resistance gene (co-8) has been reported in the common bean accession AB 136. However, when we performed PCR amplification with our tightly linked marker IAC 238, we got different amplicons in AB 136 and KRC 8. Interestingly, the susceptible cultivar Jawala produced the same amplicon as AB 136. This observation indicated that the recessive gene present in KRC 8 is different from co-8. As the gene is located far away from the Co-1 locus, we suggest naming the recessive gene co-Indb/co-19. Fine mapping of co-Indb in KRC 8 may provide new insights into the cloning and characterization of this recessive gene so that it can be incorporated into future bean improvement programs. Further, the tightly linked marker IAC 238 can be utilized in marker assisted introgression in future bean breeding programs. CONCLUSION: The novel co-Indb gene present in Himalayan landrace KRC 8, showing adult plant resistance against common bean anthracnose, is independent from all the resistance genes previously located on chromosome Pv01.

Consistent effects of independent domestication events on the plant microbiota.

The effect of plant domestication on plant-microbe interactions remains difficult to prove. In this study, we provide evidence of a domestication effect on the composition and abundance of the plant microbiota. We focused on the genus Phaseolus, which underwent four independent domestication events within two species (P. vulgaris and P. lunatus), providing multiple replicates of a process spanning thousands of years. We targeted Phaseolus seeds to identify a link between domesticated traits and bacterial community composition as Phaseolus seeds have been subject to large and consistent phenotypic changes during these independent domestication events. The seed bacterial communities of representative plant accessions from subpopulations descended from each domestication event were analyzed under controlled and field conditions. The results showed that independent domestication events led to similar seed bacterial community signatures in independently domesticated plant populations, which could be partially explained by selection for common domesticated plant phenotypes. Our results therefore provide evidence of a consistent effect of plant domestication on seed microbial community composition and abundance and offer avenues for applying knowledge of the impact of plant domestication on the plant microbiota to improve microbial applications in agriculture.

Comprehensive Proteomic Analysis of Common Bean (Phaseolus vulgaris L.) Seeds Reveal Shared and Unique Proteins Involved in Terminal Drought Stress Response in Tolerant and Sensitive Genotypes.

This study identified proteomic changes in the seeds of two tolerant (SB-DT3 and SB-DT2) and two sensitive (Merlot and Stampede) common bean genotypes in response to terminal drought stress. Differentially expressed proteins (DEPs) were abundant in the susceptible genotype compared to the tolerant line. DEPs associated with starch biosynthesis, protein-chromophore linkage, and photosynthesis were identified in both genotypes, while a few DEPs and enriched biological pathways exhibited genotype-specific differences. The tolerant genotypes uniquely showed DEPs related to sugar metabolism and plant signaling, while the sensitive genotypes displayed more DEPs involved in plant-pathogen interaction, proteasome function, and carbohydrate metabolism. DEPs linked with chaperone and signal transduction were significantly altered between both genotypes. In summary, our proteomic analysis revealed both conserved and genotype-specific DEPs that could be used as targets in selective breeding and developing drought-tolerant common bean genotypes.

Translocations and inversions: major chromosomal rearrangements during Vigna (Leguminosae) evolution.

KEY MESSAGE: Inversions and translocations are the major chromosomal rearrangements involved in Vigna subgenera evolution, being Vigna vexillata the most divergent species. Centromeric repositioning seems to be frequent within the genus. Oligonucleotide-based fluorescence in situ hybridization (Oligo-FISH) provides a powerful chromosome identification system for inferring plant chromosomal evolution. Aiming to understand macrosynteny, chromosomal diversity, and the evolution of bean species from five Vigna subgenera, we constructed cytogenetic maps for eight taxa using oligo-FISH-based chromosome identification. We used oligopainting probes from chromosomes 2 and 3 of Phaseolus vulgaris L. and two barcode probes designed from V. unguiculata (L.) Walp. genome. Additionally, we analyzed genomic blocks among the Ancestral Phaseoleae Karyotype (APK), two V. unguiculata subspecies (V. subg. Vigna), and V. angularis (Willd.) Ohwi & Ohashi (V. subg. Ceratotropis). We observed macrosynteny for chromosomes 2, 3, 4, 6, 7, 8, 9, and 10 in all investigated taxa except for V. vexillata (L.) A. Rich (V. subg. Plectrotropis), in which only chromosomes 4, 7, and 9 were unambiguously identified. Collinearity breaks involved with chromosomes 2 and 3 were revealed. We identified minor differences in the painting pattern among the subgenera, in addition to multiple intra- and interblock inversions and intrachromosomal translocations. Other rearrangements included a pericentric inversion in chromosome 4 (V. subg. Vigna), a reciprocal translocation between chromosomes 1 and 5 (V. subg. Ceratotropis), a potential deletion in chromosome 11 of V. radiata (L.) Wilczek, as well as multiple intrablock inversions and centromere repositioning via genomic blocks. Our study allowed the visualization of karyotypic patterns in each subgenus, revealing important information for understanding intrageneric karyotypic evolution, and suggesting V. vexillata as the most karyotypically divergent species.

Fine mapping of a new common bean anthracnose resistance gene (Co-18) to the proximal end of Pv10 in Indian landrace KRC-5.

KEY MESSAGE: Mapping and fine mapping of bean anthracnose resistance genes is a continuous process. We report fine mapping of anthracnose resistance gene Co-18 which is the first anthracnose gene mapped to Pv10. The discovery of resistance gene is a major gain in the bean anthracnose pathosystem research. Among the Indian common bean landraces, KRC-5 exhibit high levels of resistance to the bean anthracnose pathogen Colletotrichum lindemuthianum. To precisely map the anthracnose resistance gene, we used a Recombinant Inbred Line (F2:9 RIL) population (KRC-5 × Jawala). The inheritance test revealed that KRC-5 carries a dominant resistance gene temporarily designated as Co-18. We discovered two RAPD markers linked to Co-18 among 287 RAPD markers. These RAPD markers were eventually developed into SCARs (Sc-OPR15 and Sc-OPF6) and flank Co-18 on chromosome Pv10 at a distance of 5.3 and 4.2 cM, respectively. At 4.0-4.1 Mb on Pv10, we detected a SNP (single-nucleotide polymorphism) signal. We synthesized 58 SSRs and 83 InDels from a pool of 135 SSRs and 1134 InDels, respectively. Five SSRs, four InDels, and two SCARs were used to generate the high-density linkage map, which led to the identification of two SSRs (SSR24 and SSR36) that are tightly linked to Co-18. These two SSRs flank the Co-18 to 178 kb genomic region with 13 candidate genes including five NLR (nucleotide-binding and leucine-rich repeat) genes. The closely linked markers SSR24 and SSR36 will be used in cloning and pyramiding of the Co-18 gene with other R genes to develop durable resistant bean varieties.

Possible melatonin-induced salt stress tolerance pathway in Phaseolus vulgaris L. using transcriptomic and metabolomic analyses.

Melatonin plays important roles in multiple stress responses; however, the downstream signaling pathway and molecular mechanism remain unclear. This study aimed to elucidate the transcriptional regulation of melatonin-induced salt stress tolerance in Phaseolus vulgaris L. and identify the key downstream transcription factors of melatonin through transcriptomic and metabolomic analyses. The melatonin-induced transcriptional network of hormones, transcription factors, and functional genes was established under both control and stress conditions. Among these, eight candidate transcription factors were identified via gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, one gene related to transmembrane transport of salts (Phvul.004G177300). These genes may play a role in maintaining the cell structure and excreting sodium ions outside the cell or transporting them to the vacuoles for storage. Melatonin regulates the Phvul.009G210332 gene and metabolites C05642 (N-acetyl-N-2-formyl-5-methoxycanurine), C05643 (6-hydroxymelatonin), C05660 (5-methoxyindoleacetic acid) involved in tryptophan metabolism. The metabolites C05642 and C05643 were identified as decomposition products of tryptophan, indicating that exogenous melatonin entered the P. vulgaris tissue and was metabolized. Melatonin promotes the synthesis and metabolism of tryptophan, which is crucial to plant metabolism, growth, maintenance, and repair.

Molecular Characteristics of Bean Common Mosaic Virus Occurring in Inner Mongolia, China.

Bean common mosaic virus (BCMV) was detected on common bean (Phaseolus vulgaris) plants showing wrinkled and/or narrow leaves, curling, shrinking and chlorosis of leaves, dwarfing of plants, and mottled pods in Inner Mongolia and named BCMV-22Huhe. Its genome has a size of 10,062 bp and was deposited in GenBank under the accession number OR778613. It is closely related to BCMV-Az (GenBank accession no. KP903372, in China) in the lineage of AzBMV. A recombination event was detected for BCMV-22Huhe among the 99 BCMV isolates published in the NCBI GenBank database, showing that BCMV-CJ25 (MK069986, found in Mexico) was a potential major parent, and the minor parent is unknown. This work is the first description of the occurrence of BCMV in Inner Mongolia, China.

Genome-Wide Association Study Reveals a QTL Region for Ashy Stem Blight Resistance in PRA154 Andean Common Bean.

Ashy stem blight (ASB) caused by Macrophomina phaseolina (Tassi) Goidanich affects the common bean (Phaseolus vulgaris L.) at all growing stages. Higher levels of resistance were observed in Andean common beans, but specific resistant quantitative trait loci (QTLs) conferring resistance to this pathogen have not been reported in this gene pool. The objectives of this research were to: (i) conduct a genome-wide association study (GWAS) and QTL mapping for resistance in the Andean breeding line PRA154; and (ii) identify single nucleotide polymorphism (SNP) markers and candidate genes for ASB resistance. Phenotyping was conducted under greenhouse conditions by inoculating the 107 F6:7 recombinant inbred lines (RILs) derived from the cross between the susceptible cultivar 'Verano' and the partial-resistant breeding line PRA154 twice with the M. phaseolina isolate PRI21. Genotyping was performed with 109,040 SNPs distributed across all 11 P. vulgaris chromosomes. A novel major QTL was located between 28,761,668 and 31,263,845 bp, extending 2.5 Mbp on chromosome Pv07, and the highest significant SNP markers were Chr07_28761668_A_G, Chr07_29131720_G_A, and Chr07_31263845_C_T with the highest LOD (more than 10 in most of the cases) and R-squared values, explaining 40% of the phenotypic variance of the PRI21 isolate. The gene Phvul.007G173900 (methylcrotonyl-CoA carboxylase alpha chain and mitochondrial 3-methylcrotonyl-CoA carboxylase 1 [MCCA]) with a size of 10,891 bp, located between 29,131,591 and 29,142,481 bp on Pv07, was identified as one candidate for ASB resistance in PRA154, and it contained Chr07_29131720_G_A. The QTL and genetic marker information could be used to assist common bean breeders to develop germplasm and cultivars with ASB resistance through molecular breeding.

A novel pectin methylesterase inhibitor, PMEI3, in common bean suggests a key role of pectin methylesterification in Pseudomonas resistance.

The mechanisms underlying susceptibility to and defense against Pseudomonas syringae (Pph) of the common bean (Phaseolus vulgaris) have not yet been clarified. To investigate these, 15-day-old plants of the variety Riñón were infected with Pph and the transcriptomic changes at 2 h and 9 h post-infection were analysed. RNA-seq analysis showed an up-regulation of genes involved in defense/signaling at 2 h, most of them being down-regulated at 9 h, suggesting that Pph inhibits the transcriptomic reprogramming of the plant. This trend was also observed in the modulation of 101 cell wall-related genes. Cell wall composition changes at early stages of Pph infection were associated with homogalacturonan methylation and the formation of egg boxes. Among the cell wall genes modulated, a pectin methylesterase inhibitor 3 (PvPMEI3) gene, closely related to AtPMEI3, was detected. PvPMEI3 protein was located in the apoplast and its pectin methylesterase inhibitory activity was demonstrated. PvPMEI3 seems to be a good candidate to play a key role in Pph infection, which was supported by analysis of an Arabidopsis pmei3 mutant, which showed susceptibility to Pph, in contrast to resistant Arabidopsis Col-0 plants. These results indicate a key role of the degree of pectin methylesterification in host resistance to Pph during the first steps of the attack.

Identification of novel sources of partial and incomplete hypersensitive resistance to rust and associated genomic regions in common bean.

Common bean (Phaseolus vulgaris) is one of the legume crops most consumed worldwide and bean rust is one of the most severe foliar biotrophic fungal diseases impacting its production. In this work, we searched for new sources of rust resistance (Uromyces appendiculatus) in a representative collection of the Portuguese germplasm, known to have accessions with an admixed genetic background between Mesoamerican and Andean gene pools. We identified six accessions with incomplete hypersensitive resistance and 20 partially resistant accessions of Andean, Mesoamerican, and admixed origin. We detected 11 disease severity-associated single-nucleotide polymorphisms (SNPs) using a genome-wide association approach. Six of the associations were related to partial (incomplete non-hypersensitive) resistance and five to incomplete hypersensitive resistance, and the proportion of variance explained by each association varied from 4.7 to 25.2%. Bean rust severity values ranged from 0.2 to 49.1% and all the infection types were identified, reflecting the diversity of resistance mechanisms deployed by the Portuguese germplasm.The associations with U. appendiculatus partial resistance were located in chromosome Pv08, and with incomplete hypersensitive resistance in chromosomes Pv06, Pv07, and Pv08, suggesting an oligogenic inheritance of both types of resistance. A resolution to the gene level was achieved for eight of the associations. The candidate genes proposed included several resistance-associated enzymes, namely ß-amylase 7, acyl-CoA thioesterase, protein kinase, and aspartyl protease. Both SNPs and candidate genes here identified constitute promising genomics targets to develop functional molecular tools to support bean rust resistance precision breeding.

Genetic control of pod morphological traits and pod edibility in a common bean RIL population.

KEY MESSAGE: QTL mapping, association analysis, and colocation study with previously reported QTL revealed three main regions controlling pod morphological traits and two loci for edible pod characteristics on the common bean chromosomes Pv01 and Pv06. Bean pod phenotype is a complex characteristic defined by the combination of different traits that determine the potential use of a genotype as a snap bean. In this study, the TUM RIL population derived from a cross between 'TU' (dry) and 'Musica' (snap) was used to investigate the genetic control of pod phenotype. The character was dissected into pod morphological traits (PMTs) and edible pod characteristics (EPC). The results revealed 35 QTL for PMTs located on seven chromosomes, suggesting a strong QTL colocation on chromosomes Pv01 and Pv06. Some QTL were colocated with previously reported QTL, leading to the mapping of 15 consensus regions associated with bean PMTs. Analysis of EPC of cooked beans revealed that two major loci with epistatic effect, located on chromosomes Pv01 and Pv06, are involved in the genetic control of this trait. An association study using a subset of the Spanish Diversity Panel (snap vs. non-snap) detected 23 genomic regions, with three regions being mapped at a position similar to those of two loci identified in the TUM population. The results demonstrated the relevant roles of Pv01 and Pv06 in the modulation of bean pod phenotype. Gene ontology enrichment analysis revealed a significant overrepresentation of genes regulating the phenylpropanoid metabolic process and auxin response in regions associated with PMTs and EPC, respectively. Both biological functions converged in the lignin biosynthetic pathway, suggesting the key role of the pathway in the genetic control of bean pod phenotype.